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The probability distribution of the entropy production for the effusion of an ideal gas between 
two compartments is calculated explicitly. The fluctuation theorem is verified. The analytic results 
are in good agreement with numerical data from hard disk molecular dynamics simulations. 
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I. INTRODUCTION 



Effusion is the motion of a gas through a small pore or 
opening. The linear dimension of the opening is assumed 
to be smaller than the mean free path of the gas parti- 
cles. Consequently, the leak does not disturb the state 
of the gas. If the latter is at equilibrium, it will remain 
so and many properties of the effusing gas particles can 
be easily calculated, using basic arguments from kinetic 
theory of gases. One of the older results is the so-called 
Graham's law p), stating that the rate of effusion is in- 
versely proportional to the square root of the mass of 
the gas particles. This observation was in fact used in 
the production of the atomic bom during the Manhattan 
project for the enrichment of uranium. Graham's law is 
a direct consequence of the obvious fact that faster par- 
ticles will exit through the hole more frequently. In fact, 
the average kinetic energy of the escaping particles turns 
out to be 2fcT, to be compared to the bulk average value 
of (3/2)fcT (3 dimensional system). As a result, effusion 
will lead to a cooling of the gas, a principle that is used 
to date for refrigeration Q. In the early 1900's, Knudsen 
investigated other aspects of effusion 0]. He discovered 
the cosine law, corresponding to the fact that the parti- 
cles exit with an isotropic angular distribution, a prop- 
erty which was subsequently using for coating spherical 
bulbs, and which is also related to the Kirchoff 's laws of 
radiation and to proper boundary condition for particles 
that are thermally reemitted from an absorbing wall |4| . 

An interesting situation occurs when the small hole 
forms the contact between two containers (referred to as 
Knudsen cells) with gases which are separately at equi- 
librium, but not at equilibrium with each other, cf. Fig. 
^ While local equilibrium (inside each container) is pre- 
served, the existing density and temperature gradients, 
which can be arbitrarily large, will lead to both particle 
and energy fluxes. A somewhat surprising observation, 
already noted by Knudsen, is that for paVTa = PbVTb 
an energy flux is present in the absence of a correspond- 



ing particle flux, cf. Eq. (|5U[) (this peculiar property also 
appears in the so-called Jepsen gas |5|). Similarly, a par- 
ticle flux is present with no corresponding energy flux 
when p a Ta 12 = p B T B 3/2 , cf. Eq. g5J. When oper- 
ating in the regime of weak gradients, more precisely in 
the framework of linear irreversible thermodynammics, 
the properly defined thermodynamic forces and fluxes 
are related to each other by a symmetric Onsager ma- 
trix labUs!- 

The above properties refer to average values of the 
fluxes. It is however clear that one can easily calculate 
fluctuations and correlation functions. For example, the 
typical appearance of long range spatial correlations in 
nonequilibrium systems was illustrated in detail for an 
array of Knudsen cells [1 [13 [H G3 . 

In the present paper, we go one step further and cal- 
culate the full joint probability distribution for the par- 
ticle and energy flux. This result is particularly relevant 
because of the recent discover y o f various work and fluc- 
tuation theorems E3, EM [IH E3, EH E3, E3, H3, 

I23I |24| , which result from a time reversal symmetry of 
the underlying dynamics. In the present case, it implies 
that the probability density Pt(AS) for the entropy AS, 
produced by effusion during a time interval of duration 
t, satisfies the following relation: 



ft (AS) 
Pt(-AS) 



= e AS/k . 



This result implies: 



-AS/k\ 



(1) 



(2) 
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and consequently (invoking Jensen's inequality) (AS) > 
0, in agreement with the second law. 

The purpose of this paper is to present an explicit and 
detailed analysis of effusion in the context of various fluc- 
tuation theorems. Their verification is most easily done 
at the level of the cumulant generating function. The ef- 
fusion problem moreover has the pleasant feature that the 
marginal distribution functions for particle and energy 
transport can also be calculated explicitly. They obey 
separate fluctuation theorems when the thermodynamic 
force for the other transport process vanishes. Their large 
deviation functions can be calculated explicitly and are 



2 




P A ' T A Pb> t b 



FIG. 1: The system under consideration: two ideal gases, 
in equilibrium at their respective temperatures and densities, 
can exchange energy and particles through a small hole in the 
common adiabatic wall. 



related by Legendre transformation to the cumulant gen- 
erating function. Finally, we report the explicit results 
for the two-dimensional version, allowing a detailed com- 
parison with numerical data from hard disk molecular 
dynamics simulations. 



II. FLUCTUATION THEOREM FOR EFFUSION 

Consider two large reservoirs, A and B, separated by 
a common adiabatic wall, each containing an ideal gas 
at equilibrium, with uniform density pi and Maxwellian 
velocity distribution 4>i(v) at temperature Tj, i S {A, B} 
(cf. Fig.P): 



(3) 



At time t = , a small hole with surface area a is opened, 
allowing an exchange during a time interval of length t 1 
of particles and energy between the two compartments. 
The dimensions of the hole are assumed small compared 
to the mean free path of the gas particles, and the reser- 
voirs sufficiently large so that the equilibrium state in 
both parts is not affected by the exchange. In particular, 
the temperatures and densities in both reservoirs remain 
constant j^. Hence, the change in entropy in the to- 
tal system, upon a total transfer of energy AU and of 
particles A7V from A to B during t, is given by 



AS 



= AS A 
1 

~T~A 
= AuAU 



AS E 



AU- 



PA 



AN 



■A N AN. 



(4) 



We introduced, in accordance with the definitions from 
irreversible thermodynamics, the following thermody- 
namic forces (affinities) for energy and particle flow, re- 



spectively: 
A v 



1 

Tb~ 



1 

T~A 



A 



N 



HA PB 

77 TTT - = fclOg 



PA 
PB 



Tb 
T a 



(5) 



In the last equality, we used the expression for the chem- 
ical potential of an ideal gas. Note that the thermo- 
dynamic force An diverges when one of the (infinitely 
large!) reservoirs is empty, so that free effusion into un- 
limited space, implying an infinitely large entropy pro- 
duction, corresponds to a singular limit. 

Clearly, as single particle crossings from both sides 
contribute to the values of AU and AN in the course 
of time, these quantities are fluctuating, and hence so is 
the entropy production AS. As indicated in the intro- 
duction, the resulting probability density for the entropy 
production should obey the fluctuation theorem, Eq. 
Clearly, the contributions to AS from any two equal but 
non-overlapping time intervals are independent and iden- 
tically distributed random variables. In other words, AS 
is a stochastic process with independent increments. As a 
result its cumulant generating function has the following 
form: 



(6) 



The fluctuation theorem Eq. Q implies the following 
symmetry relation for the function /i(A) (to which we 
will also refer, for brevity, as the cumulant generating 
function) : 



p(X) = Mfc- 1 - A). 



(7) 



We will in fact prove for the system under consideration 
a more detailed fluctuation theorem, expressed in terms 
of the joint probability density, namely: 



P t (AU, AN) 
P t (-AU,-AN) 



a AS/k 



(8) 



The increments of AU and AiV are again independent 
so that the corresponding cumulant generating function 
has the form: 



i e -{\uAU+\ N AN)\ _ e -tM(A[/,A N ) 



(9) 



In terms of the function p{Xjj , Xn), the detailed fluctua- 
tion theorem reads: 



n(\u, A w ) = p,{Au/k - Xu,AN/k - Ajv). 



(10) 



Since AS" = AuAU + A N AN, one has fi(X) = 
p(XAu, XAn), so that the detailed fluctuation theorem 
implies the "normal" fluctuation theorem. 

Eq. (JSJ also implies fluctuation theorems for energy 
and particle transport separately, when the correspond- 
ing thermodynamic force for the other process vanishes. 
Indeed, in this case only one of the variables, AU or AN, 



3 



appears in AS. Bringing the denominator from lhs to rhs 
in Eq. (JSJ , subsequent integration over the other variable 
implies: 



Vt(AU) 
Vt(-AU) 



= e AS/k when A N = 0, (11) 



and 



VtjAN) 
Vt(-AN) 



3 AS/k when Au = 0. (12) 



We now proceed to a direct verification of the detailed 
fluctuation theorem Eq. I|10|l and its implications Eqs. 
()7I11I12(I by an explicit evaluation of the generating func- 
tion ii(\u,X n ). 



III. MASTER EQUATION AND CUMULANT 
GENERATING FUNCTION 



During a small time interval dt, the contributions to 
quantities AU and AN result from individual particle 
transport across the hole. Following basic arguments 
from kinetic theory of gases (for more details, see Ap- 
pendix^), the corresponding probabilities per unit time, 
Ta^b(E) and Tb^a(E), to observe a particle with ki- 
netic energy ^mv 2 — E crossing the hole from A — > B 
and B — > A respectively, are given by: 



ap A E 



yJ2nmkT A kT A 
apB E 



e kT A 



^2wmkT B kT B 



g K1 B 



(13) 



The probability density P t (AU, AN) thus obeys the fol- 
lowing Master equation: 

J^; -P* (At/, AN) = J T A ^B(E)P t (AU-E,AN-l)dE 

+ / T B ^A(E)P t (AU + E,AN + l)dE 
Jo 

/>oo 

-P t (AU,AN) [T A ^B{E)+T B ^ A {E)}dE. (14) 
Jo 

Note that the integral operators are of the convolution 
type, in agreement with the fact that the processes AU 
and AN have independent increments. Hence, the ex- 
act solution of Eq. Q14JI. subject to the initial condi- 
tion Po(AU, AN) — S(AU)San,o, is obtained by Fourier 
transfom, i.e., by switching to the cumulant generating 
function /j,(\u, Xn). Since the following integrals can be 
performed explicitly: 



a paVWa 



/ /, .,;</•;< A ' 1 di: v 

Jo V27TTO (1 + kT A Xu) 

[°°T B ^ A (E)e^ E dE = a PbVWE 
Jo 



V2^ (1 - kT B \uY 



, (15) 



we immediately obtain: 

p(Xu,X N ) = J y== ( /'.i \ T \ 



+Pb\/Tb 



1 - 



(1 + kT A XuY 

Xn 



(1 - kT B XuY 



(16) 



Eq. (|16|l is a central result of this paper. One easily veri- 
fies, using the explicit expressions for the thermodynamic 
forces given in Eq. J^J), that this expression indeed verifies 
the detailed fluctuation theorem Eq. (|1(J|) . 

Note also that p,(Xu, Xn) can be written as the sum of 
two contributions: 



p(Xu,Xn) — Pa(^u 7 Ajv) + Hb{—Xu, —Xn), 

with 

P-a(Xu, Xn) 



apAVkTA 



V2r 



1 



(1 + kTAXu)- 2 



(17) 



(18) 



and the similar expression for A changed into B. This 
additivity property derives from the statistical indepen- 
dence of the fluxes from A — > B and B — > A. 



IV. CUMULANTS 

The joined cumulant Ky of power i in energy flux and 
j in particle flux appears as a coefficient of the Taylor 
expansion of the cumulant generating function, namely: 



i,j=0 



ilj\ 



The explicit expressions for the cumulants of energy and 
particle flux can thus easily be obtained from Eq. (TBJ. 
We mention the following results. The first order cumu- 
lants read: 



K W = (AU) = ^2 ( Pa tT - PB Tf 2 ) ; (20) 
«oa = (AN) = t -^l(pAT^-p B T^). (21) 



The second order cumulants are 



n 2Q = (SAU 2 ) = ^6 (paTT + p B T 5 J 2 ) : ( 22 ! 



tak 3 / 2 



til 



k Q2 = (SAN 



2 [paT Z a 12 + PB T 3 B /2 ) ; 



(23) 



; tak 1 / 2 / 1/2 1/2 \ 

= -7= \PaT a ! +PbT b ' ) . (24) 
V27rm v 7 



Similar expressions can be obtained for the higher order 
cumulants. Note that all cumulants are proportional to 
time, the specifying feature of processes with indepen- 
dent increments. 
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V. MARGINAL DISTRIBUTIONS 

The explicit calculation of P t (AU,AN) is quite in- 
volved, but it is relatively easy to obtain analytic ex- 
pressions for the marginal distributions: 

/+oo 
P t (AU,AN)dAU; (25) 
- CO 
+00 

Vt(AU) = Yl Pt(&U,AN). (26) 

AJV=-oo 

The corresponding cumulant generating functions of 
these marginal distributions are obtained by setting 
A at = and \u = respectively in Eq. fTH|l . 

The generating function for particle transport is thus 
found te be: 

,(0,A)^( M ^[l-e1 
V Zirm v 

+ p B Vf^[l-e x ]y (27) 
This result is identical to that for a random walk taking 



steps to the right and left with probabilities per unit time: 

r A = apA\ZkTA/2-Km; (28) 
y/kT B /2irm, (29) 

respectively. This relation is intuitively clear since the 
passages of particles from reservoir A to B and vice- versa 
induce a random walk on the variable AN. A similar 
expression for the cumulant generating function, but with 
re = 0, appears in the problem of shot noise [25t |26|. 
The probability distribution P t {AN) is, as is well known 
from the random walk literature, expressed in terms of 
the modified Bessel function. Since we have: 

e -T/*(o,A) = J- e- XAN V t (AN), (30) 

AN 

and using the generating function of the modified Bessel 
functions |27| : 

e M e - A +^) = ^ e -Ar %(z)j (31) 

n 

one finds: 



V t {AN) 



• ^PA^kTA+PB 




(paPbVTaTb) 1 



/2 



(32) 



When Ta — Tb, i.e., when the thermodynamic force for 
the energy flow vanishes, Au = 0, we can verify explicitly 
the fluctuation theorem Eq. JT2J: 



VtjAN) 
V t (-AN) 



PA 
PB 



AN 



AN 



^AS/k 



(33) 



The calculation of the marginal probability distribu- 
tion for the energy transport is slightly more involved. 
We will invoke the independence of fluxes from each reser- 
voir to the other one, as already alluded to when dis- 
cussing the additivity property of the generating func- 
tion. The total amount of energy AU transfered during 
a certain time interval can be written as the difference of 
two independent contributions AUa and AUb, being the 
energy coming from compartment A and B respectively: 



AU = AU A - AU E 



(34) 



It follows that Vt{AU) can be expressed as a convolution 
integral: 



Vt(AU) 



Vf{E)Vt (E- AU)dE, 



(35) 



with Pf(AU) and V^(AU) being the probability dis- 
tributions to have an energy transfer AU as a result of 
particles solely going from A — > B and B — > A, respec- 
tively. From: 



-*Ma(A[/,0) _ 



+00 



-XuAUqjA 



T t A {AU)dAU, 



(36) 



the calculation of T'f(AU) involves an inverse Fourier 
transform, which is reproduced in Appendix iBl The final 
result reads: 



rf(Au) 



-top A \ 



kT A 



5{AU) + 9(AU) 



tap a AU 



yJ2irmkT A kT A 



■ kT A 



{}, 



tapAVkTA 
4\/27rm 



'AU' 



(37) 
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The expression for V^{AU) is obtained by replacing A 

3/2 3/2 

by B in the above equation. When T A / pa = TJ / p B , 
i.e., the thermodynamic force for the particle flow van- 
ishes, An = 0, the resulting probability density Vt{AU) 
satisfies the fluctuation theorem Eq. H12|) . This follows 
directly by verifying from Eq. (|37|l that Vf-(AU) and 
Tt{MJ) obey the relation (VE 1 ): 



V^{E)Vf {E — AU) 
V t A (E-AU)V t B (E) 



^AS/k 



(38) 



is a convex function, independent of t, and known as the 
large deviation function. This function can also be ob- 
tained as the Legendre transform of the generating func- 
tion /i(0, A): 



<p{n) ~ sup{/i(0, A) — An}. 

A 



(41) 



VI. LARGE DEVIATION FUNCTION 



Using following asymptotic expansion of the Bessel func- 
tions for large orders p7|: 



The particle transport AN is expected to grow propor- 
tional to time t. The asymptotic behaviour of Vt(AN) 
for large time t — > 00 is indeed described in terms of the 
variable n = AN /t, namely: 



Vt(AN) 



-tip(r. 



(39) 




r^j (i + z 2 )i/4 



More precisely, the following limit: 



1 



ip(n) = — lim - hiVtint), 

t — >oo t 



and substituting v = nt and z = we obtain 

(40) from Eq. lj3l|) : 



V t {nt) 



-t(r A +r B ) p nt In ^Jr A /r B . 



t ( V ^ r Af B +n 2 +n In 



(4r A r B + n 2 ) 1 / 4 



(43) 



This result leads to the following expression for <p(n): 



ip(n) = ta + r B - V -irATB + n 2 



- nh \ 2^ J' (44) 

which is of course the well-know large deviation function 
for a (biased) random walk. The same result is obtained 
using the Legendre transform of the generating function 
p(0,X). A sketch of both u(0, A) and <p{n) is given in 
Fig. (|2J) - The function (p(n) is positive everywhere, and 
has a single zero at the most probable value n = n = 
ta - r B - 

Similarly, the asymptotic properties for t — > 00 of the 
energy transport are obtained by focusing on the scaled 
variable u = AU/t. The large deviation function of 
Vf(AU), 



ip A {u) 



lim -]nVf{ut), 

t^oo t 



(45) 



is again obtained either by Legendre transform of 
Pa(^, 0), cf. Eq. (Hi, or directly from Eq. (|37|) in which 
case the asymptotic behaviour of o^2({}, {3/2, 2}, z) is 




FIG. 2: Sketch of the generating function /i(0, A) and its Leg- 
endre transform, the large deviation function f(n) for the case 



needed: 

F 2 ({},{3/2,2},z) 
The final result reads: 



4V37TZ 5 / 6 



for 



(46) 



<Pa(u) 



q / 2 \ 1/3 

3 / apAU N 



kT A VkTA \4v / 2^ 



(47) 

A sketch of the generating function //^(A, 0) and (Pa(u) 
is given in Fig.@. (Pa(u) has a unique zero at the most 
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<g>) 



X 2kT A r A 



FIG. 3: Sketch of the generating function Pa(X) and the large 
deviation function ipa(u). 



where we define the macroscopic fluxes Ju and J7jv corre- 
sponding to energy and particle transport, cfr. Eqs. H20|) - 

jaj: 



d /A „, a2k 3 / 2 
Ju = —(AC/) = -—= 

Ut V27TTO 



Jiv = ^t(A^) = 

dr 



(paT^ 2 - p s ^ /2 ) ; (49) 

^ (p.tT-PbTT)- (50) 



probable value u — u — 2a p A(kJ a)^ /^/2nm — 2VTaTa- 
This result is in agreement with the fact, cf. introduction, 
that particles carry on average 2kTA energy. Note that 
T'f-(AU) decays exponentially for large AC/, explaining 

the divergence of pa{X, 0) for A ^ — 1/(/cTa), and the 
convergence of <pa{u) to a linear function for u — > 00. 



VII. RECIPROCITY RELATIONS 

Averaging Eq. (@J and taking the time derivative gives 
the average entropy production: 



— (AS)=J U A U + J N A N , 
at 



(48) 



In general, these fluxes are functions of the affinities 
(-4(7) ■Aw)- Close to equilibrium, i.e. for small tempera- 
ture and density differences AT and Ap, we can write: 

Ta ~T — AT/2 ; p A = p - Ap/2- 
T B =T + AT/2 ; p B = p + Ap/2, (51) 

so that the thermodynamic forces become: 



Av*-§ and A„ (52) 



A Taylor expansion of Ju and Jjv leads to the following 
results: 



Ju = -^L=2k 1 / 2 pT 3 / 2 (ikTAu + A N + ^kT 3 A v + ^T 2 AfjA N + ^kT 5 A v + ^T 4 ^^ +...), (53) 



and 



J N = -^2fc 1 / 2pT 3/2 ( Au + 1 

V27rm V 2Tfc 64 



1 > , IT. 
-Av - — T .Ay — — — A v An 



1024 



T A A% 



5 T 3 
4096 T 



U-^N 



(54) 



r 



From these expressions, it is clear that the Onsager sym- 
metry [HOI 13 is satisfied: 



dJu 



dA 



N 



dJi 



N 



dA, 



(55) 



where the derivatives are evaluated at equilibrium 
(Au,An) — (0,0). This relation is in fact also a di- 
rect consequence of th e sy mmetry relation of p(\u, Xn), 
cfr. Ea. (fTHft (11 I2I l29l l3fj|. To show this explicitly 
it is convenient to include the dependence of p on the 
thermodynamic forces, i.e. we write p(X, A) (for ease of 
notation, we use A = {Xu, Ajv} and A = {Au,An})- We 



then find 
dJu 



8A N 



d 2 p{X,A) 



d 2 p(A/k-X,A) 



OAnOXu 
1 d 2 p(X,0) 



k OXuOXn 
The last term is again 

dJu 



8Xu8An 
d 2 p{-X,A) 



dA N 



dXudAN 
dJu/dAN\ , so that: 
1 d 2 p(X,Q) 



2k 8XuOXn 



(56) 
(57) 

(58) 



which leads to the Onsager symmetry relation. 

VIII. 2D RESULTS 

All the above results have been derived for three spatial 
dimensions. For comparison with computer simulations, 
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we reproduce the corresponding results for two spatial 
dimensions. We note the following two main modifica- 
tions: 



i) The thermodynamic force corresponding to 
particle flux has the following form: 



the 



A 



N 



k log 



PbTa 



ii) The transition rates Ta^b(E) and T b ~>a{E) 
are to be calculated using the two-dimensional 
Maxwellian velocity distribution: 



2nkT, 



(59) namely: 



(60) 



T A ^ B {E) = 



and 



OO p77 



vdvdtp pa<I>a(v) u v sin tp 5(E — —mv 2 ) = X / ' i '- \ 



v=0 Jtp=0 



-K\JmkTA \kTA 



e kT A ; 



(61) 



OO />27T 



vdvdtp pb4>b{v) a (— v sin tp) 5(E — -mv 2 ) 



V2ap B ( E 



V—0 J (p =7T 



(62) 



r 



They display a prefactor proportional to s/E, hence tran- 
sitions with E < kT gain more weigth as compared to the 
3D situation. The calculation of the generating function 
proceeds in exactly the same way. The final result differs 
from the 3D case only by the 3/2 power law exponent: 



a\fk 




At(A[/, Aj\r) 



The fluctuation relation Eq. (|10|l is again verified. The 
first and second order cumulants read: 



«10 



(AC/) = 
koi = (AAT) = 



3 iafc 3 / 2 



2 v27rm 
tay/k 
\T2jrm 



(paT 



3/2 „3/2 

A - PbT b 



(paT a /2 - PbT 



1/2 
B 



k 2 q = {SAU ) = - — j== ypAT^ 



3 tak 3 '' 



2 y/2r. 



n 02 = (SAN 2 ) = 



4 V27rm 

[paTT 

tasfk 
\f2wm 



5/2 
A 



PbT b > 



PbT b 



PaT 



1/2 



,7,1/2 

PbT b 



(64) 
(65) 
(66) 
(67) 
(68) 



(Nb) disks are randomly placed in reservoir A (B). Their 
velocities are randomly sampled from the Maxwellian ve- 
locity distribution at the corresponding temperature and 
slightly rescaled in order to make the center-of-mass ve- 
locity zero in both reservoirs, and the total energy equal 
to NAkTA and NskTs in reservoir A and B respectively. 
Then, the system is allowed to relax for a certain period, 
after which the hole is opened and AU and AN are mea- 
sured. Averages are taken over 1 000 000 realizations. 
The duration time r of the simulation is measured in 
units of the average time between two particles going 
from A to B, and is related to t via: 



tapA\/kTA 



V2t 



(69) 



A comparison with the above given theoretical results for 
an ideal gas is shown in Figs. 14151617151 andlSl and Tables^ 
andllTI Agreement is extremely good. In particular, the 
fluctuation theorem is verified, cf. Eq. ^ I n Table |U 
and [HJ we have also included the numerically obtained 
value of the average (exp— AS/k). Note the deviations 
for increasing r from the exact theoretical value 1 due to 
unsufflcient sampling of the negative tails of Pt(AS). 



X. DISCUSSION 



IX. 2D MOLECULAR DYNAMICS 

Extensive molecular dynamics simulations were caried 
out for the 2D system. Two equally large reservoirs of 
total size 5000 x 200 (reflecting boundary conditions) and 
connected through a hole of diameter a = 5, contain di- 
lute gases of hard disks of diameter d = 1 and mass m = 1 . 
Initially, the hole between the reservoirs is closed and Na 



The results presented here, the fluctuation theorem for 
the effusion of an ideal gas, have a certain academic and 
didactic appeal. The underlying physics is easy to grasp, 
while the exact and detailed analytic solution completes 
the understanding. The results are not restricted to the 
regime of linear response since the gradients in density 
or temperature between the reservoirs can be arbitrarily 
large. Yet, a few words of caution are in place. First, the 
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TABLE I: Comparison between molecular dynamics simulation and theory (pa/pb = 2.0 and Ta/Tb = 2.0). 



r 






K01 




K-20 




Kll 




«02 








sim. 


theory 


sim. 


theory 


sim. 


theory 


sim. 


theory 


sim. 


theory 


sim. 


0.1 


0.122 


0.124 


0.062 


0.065 


0.407 


0.408 


0.172 


0.177 


0.127 


0.135 


0.997 


1.0 


1.280 


1.235 


0.657 


0.647 


4.214 


4.082 


1.802 


1.765 


1.358 


1.354 


0.955 


4.0 


5.178 


4.939 


2.667 


2.586 


16.904 


16.326 


7.247 


7.061 


5.456 


5.414 


0.773 


8.0 


10.393 


9.879 


5.378 


5.172 


33.619 


32.652 


14.463 


14.121 


10.893 


10.828 


0.707 



TABLE II: Comparison between molecular dynamics simulation and theory (pa/pb = 0.5 and Ta/Tb = 4.0). 
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(e-AS/kj 




sim. 


theory 


sim. 


theory 


sim. 


theory 


sim. 


theory 


sim. 


theory 


sim. 


0.1 


0.113 


0.113 


0.004 


0.0 


0.191 


0.398 


0.184 


0.188 


0.395 


0.2 


0.988 


1.0 


1.138 


1.125 


0.015 


0.0 


3.984 


3.984 


1.896 


1.875 


1.998 


2.0 


0.910 


3.0 


3.375 


3.375 


0.040 


0.0 


11.755 


11.953 


5.646 


5.625 


5.981 


6.0 


0.515 


6.0 


6.679 


6.750 


0.083 


0.0 


23.145 


23.906 


11.226 


11.250 


11.920 


12.0 


0.214 



P/AU) 




-5 5 10 15 AN 



FIG. 4: (Color online) Comparison between theory and sim- 
ulation with pa/pb = 2.0 and T A /T B = 2.0 (N A = 2N B = 
2000 and T A = 2T B = 1) for different time intervals, a) 
Plot of the marginal distribution V T (AU) (a Dirac Delta con- 
tribution at the origin is omitted), b) Plot of the marginal 
distribution V T ( A N). Dots: molecular dynamics simulations; 
full line: theoretical results (in b) they serve as a guide to the 
eye). 



Log(PjAU,AN) / P x (-AU,-AN) ) 




FIG. 5: 2D Plot of log(P T (A[/, AN)/P r (-AU, -AN)) for 
Pa/pb = 2.0, Ta/Tb = 2.0 (N A = 2N B = 2000 and T A = 
2T B = 1) and r = 1.0. 



much larger than the mean free path |3l|. Second, either 
due to the smallness of the hole or due to the ideality of 
the gas, we did not need to discuss the energy required to 
open and close the hole connecting the reservoirs. The 
size and impact of this contribution on the fluctuation 
theorem will obviously dep end on the system under con- 
sideration (33, l33l 1.3 lL 13a , l36| . Finally, fluctuation and 
work theorems can also be formulated for time-dependent 
situations, rather than the stationary state that was dis- 
cussed here. The trademark of these theorems is that one 
has to distinguish probability distributions (for entropy 
or work) for direct and time reversed schedules. These 
probabilities are typically not identical and in fact can 
describe quite different physical situations of the system, 
see for example [37]. 



present example is only a limited illustration of the fluc- 
tuation theorem. By considering effusion through a small 
hole, the system is at all times in local equilibrium. In 
this case, the fluctuation theorem can in fact be derived 
using equilibrium concepts only. However the fluctua- 
tion theorem remains valid away from local equilibrium, 
for example when the opening between the reservoirs is 



APPENDIX A: TRANSITION RATES 

The calculation of the transition rates as given in 
Eq. H13|) i is similar to that of the pressure in textbooks 
on kinetic theory. We focus here on Ta^b(E), since the 
result for Tb^a {E) is obtained from it by replacing A — > 
B. In the following, we write the velocity in spherical 
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FIG. 6: (Color online) Plot of P T (AS)/P T (-AS) for p A /p B = 
2.0 and T A /T B = 2.0 (N A = 2N B = 2000 and T A = 2T B = 1) 
for different time intervals: a) r = O.f, b) r = 1.0, c) r = 4.0, 
d) t — 8.0. The dashed line shows the theoretical result 
exp{-AS/k). 



P/AU) 
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FIG. 7: (Color online) Comparison between theory and sim- 
ulation with p A /p B = 0.5 and T A /T B = 4.0 {2N A = N B — 
2000 and T A = 4T B = 1) for different time intervals, a) Plot 
of the marginal distribution P T (AU) (a Dirac Delta contri- 
bution at the origin is ommitted). b) Plot of the marginal 
distribution V T (AN). Dots: molecular dynamics simula- 
tions; full line: theoretical results (in b) they serve as a 
guide to the eye). Note that the average particle flux is zero 
(Pax/Ta = p B \/T B ). 




FIG. 8: 2D Plot of log(P r (At/, AN)/P T (-AU, -AN)) for 
p A /p B = 0.5, Ta/Tb = 4.0 (27Va = iVs = 2000 and T A = 
4T S = 1) and r = 1.0. 
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FIG. 9: (Color online) Plot of P t (AS)/Pt(-AS) for pa/pb = 
0.5 and T A /T B = 4.0 (2N A = N B = 2000 and T A = 4T S = 1) 
for different time intervals: a) r = 0.1, b) r = 1.0, c) r = 3.0, 
d) r = 6.0. The dashed line shows the theoretical result 
exp(-A5 , /fc). 



coordinates, i.e. v — (v sin^ cos ip, v sm9 sin tp, v cos 0), 
with the z-axis perpendicular to the common wall be- 
tween the two reservoirs, and pointing from A to B, see 

Fig. uni 

TA^B(E)dEdt is the probability to observe a particle 
with kinetic energy in the range ]E,E + cLE[, crossing 
the hole from i to B during a time interval d£. Two 
requirements must be fullfilled in order to observe this 
crossing. First, the kinetic energy of the particle must 
be in the specified energy range, which means that its 
velocity is in the range ]v, v + dv[ with: 



v = y / 2E/m and dv = dE/V2Em. 



(Al) 



Second, the particle must be able to reach the hole dur- 
ing the time interval dt. A particle with velocity v cos 9 
in the z-direction must therefore be located inside the 
volume av cos Odt (cfr. Fig. llOfl . Since the z-component 
of the velocity must be positive if the particle is to move 
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from A to B; 9 has to be between and tt/2. cording to the Maxwellian velocity distribution </>a(v), 

Adding all contributions from the different directions 9 leads to the following result: 
and if and noting that the velocities are distributed ac- 



T A ^ B (E)dEdt 



it/2 r2-n 

I v 2 sin 9dvd9dtp pa4>a{v) ov cos 9dt 

=0 Jtp=0 
tt/2 ,2tt 2E 

=0 Jip=0 m 
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y/2TrmkT A kT A 
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2irkT A 
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e kT A <j\ — cosftdi 
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(A2) 




i; cos(9) dt 



FIG. 10: Sketch of the gas near the hole. A particle with 
velocity v in must be located inside the tilted cilinder with 
volume av cos 6 At, in order to reach the hole during the time 
interval At. 



APPENDIX B: CALCULATION OF THE 
MARGINAL DISTRIBUTION 

Our starting point in the calculation of the marginal 
probability distribution V^(AU) is the generating func- 
tion /xa(A,0), cf. Eq. (jl7JI . For ease of notation, we use 
the following dimcnsionlcss quantities: 



tap A VkT A Pa T a 

,- r ■■ j= iP= -X=7^, ( B1 ) 

kJ-A V27rm Pb 1b 

which are respectively the energy measured in terms of 
the thermal energy in part A, the average time between 
two particles going from A to B, the ratio of the densities 
and the ratio of the temperatures. In these new variables, 
the generating function becomes: 



AU(A,0) = 1 



1 



(1 + A) 5 



(B2) 



Taking into account the Jacobian of the transformation 
form AU -> X, and defining Vf(AU)dAU = f T (X)dX 
we obtain: 



-XX 



f T (X)dX. 



(B3) 



By inverse Fourier transform (setting A = iq), we find: 



fr(X) 



1 

27 



lqX dq. 



(B4) 



Under the integral, we add and substract the function 
e~ T e tq . The purpose for doing is twofold. First, it 
allows to single out the Dirac delta at the origin: 



um 



e {i+i q y 



1 ) J qX dq 



(B5) 

Second, the remaining integral can be converted to a con- 
tour integral using Jordan's Lemma by noting that the 
complex function exp n-^p — 1 goes to zero for large 
values of \z\. Then, for X > (X < 0), we can close the 
contour in the upper (lower) half plane. Since the inte- 
grand has only one (essential) singularity at z = i, i.e. 
in the upper half plane, and using the residue theorem it 
follows that f T (X < 0) = 0, a result that was expected 
since a particle going from A to B always takes a positive 
amount of energy with it (its kinetic energy). Using the 
Heaviside function 9{x) we conclude: 



f T (X)=e- T 5(X) + 9(X) — 



= e- T S(X) + 9(X)e- T iRes(i). 



1 



3 A" 



dz 



(B6) 



The residue of the integrandum is easily obtained by 
making a Laurent expansion around z = i: 




r) n (iX)( 2 "+J) 
(2n + j)\ 



(z-iy. (B7) 



And so: 



f T (X) = e- r S(X) + 9(X)e-^ V - 

' " 71 



n _X"(2n-l) 



^ n\ (2n- 1)! 

71=1 ^ ' 



(B8) 
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Finally, this last summation can be done explicitly in The final expression for Vf(AU) is given in Eq. (|37|l . 
terms of the generalized hypergeometric function n-Fo|39j|: 

f T (X) = e- T 5(X) 

+ 9(X)e-^rX F 2 ({}, d, 2}, . (B9) 
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